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Abstract 

The (^-dressing scheme based on local nonlinear vector 9-problcm 
is developed. It is applicable to multidimensional nonlinear equations 
for vector fields, and, after Hamiltonian reduction, to heavenly equa- 
tion. Hamiltonian reduction is described explicitely in terms of the 
<9-data. An analogue of Hirota bilinear identity for heavenly equation 
hierarchy is introduced, r-function for the hierarchy is defined. Addi- 
tion formulae (generating equations) for the r-function are found. It is 
demonstrated that r-function for heavenly equation hierarchy is given 
by the action for 9-problem evaluated on the solution of this problem. 



1 Introduction 

Dispersionless integrable equations represent themselves an important class 
of nonlinear PDEs with various applications in physics. Recently, consider- 
able interest has been paid to hydrodynamic type equations which arise in 
the general relativity, Yang-Mills theory and theory of Einstein- Weyl spaces 
(see e.g. PQ-[H]). Heavenly equations [Q have been studied particularly in- 
tensively. Different techniques have been used to analyze these equations 
and their properties. 

In the present paper we apply the 9-dressing method to such equations. 
The 9-dressing method has been proposed in ^U] and it is applicable to 
various classes of PDEs (see e.g. I12j). During the last five years it was 
demonstrated that the quasiclassical 9-dressing method is an effective tool 
to study dispersionless integrable equations (see e.g. Here we will 

show that the 9-dressing method based on local nonlinear vector 9-problem 
effectively works for class of equations for vector fields including heavenly 
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equations. Hamiltonian 3-problem generating the heavenly equations is of 
particular interest. 

2 Basic dressing scheme. Integrable equations for 
vector fields. 

We start with nonlinear local vector d-problem 

dS l = W i (z,z;S\...,S N ), i = l,...,N, (1) 

where z € C, bar means complex conjugation, 8 = J=, and functions W l 
(3-data) are smooth functions of S l . Similar to standard 9-dressing method 
[TO] , we suggest that 3-data W l (z, z; S 1 , . . . , S ) are equal to zero in some 
domain (or set of domains) G of the complex plane, thus <9-data are localized 
on a closed subset C\G. This subset may consist of a number of curves, 
thus giving rise to vector nonlocal Riemann problem. In the following we 
will usually suggest that C \ G is a unit disc D. 

In the simplest case N = 1 equation Q defines the so-called generalized 
analytic functions (see e.g. [SOD- Various classes of elliptic systems of the 
type © are also studied quite well (see e.g. |2*T] . [22] ) . 

The (^-dressing method starts with specification of analytic properties 
of solutions of the system Namely, the problem is to find the vector 
function S = (S l , . . . , S ) of the form S = So + S, satisfying Q in the unit 
disc, with S analytic outside the unit disc and decreasing at infinity, and So 
analytic in the unit disc, for arbitrary So (So plays the role of the boundary 
condition). The function S in this case may be considered as a functional 
of So- Introducing a parameterization of So, 

oo 

Sl = Y,tz n -\ (2) 

n=l 

we may also consider S as a function of iV infinite sets of times (t 1 , . . . , t^). 
First derivatives of S over times satisfy a linear system 

Bf = J2w^z,z;S\...,S N )f, (3) 

3 

where W*- = The linear space of functions, defined by this equation, 

admits multiplication by arbitrary scalar function of times and by arbitrary 
scalar functions of z, analytic inside the unit disc. 
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Solutions of equation (j3J) which are bounded on the whole complex plane 
and decrease at infinity play a special role. In the case N = 1 there is 
an analogue of the Liouville theorem (see e.g. Theorem 3.11 from |2Uj ) 
which states that a solution of (J3J) continuous and bounded on the whole 
complex plane and vanishing at some fixed point (in particular, z = oo) is 
identically zero. For systems (J5J with N > 1 this is in general not true (see 
e.g. |21|.|22j). In our approach we assume the validity of an analogue of 
Liouville theorem. It is obviously valid for the systems © with small W*-. 

However, below we will demonstrate that for some special reductions of 
the problem (^Q) it is possible to prove an analogue of Liouville theorem for 
systems (JHJ for arbitrary 3-data satisfying the reduction (given the existence 
of the basis of solutions) . 

Using the properties of the linear space of solutions of the system (J3J) and 
the standard ideology of the 5-dressing method , it is possible to 

demonstrate that the function S satisfies an (infinite) set of linear equations. 
Compatibility condition for a pair of such equations defines a closed system 
of multidimensional nonlinear PDEs for the coefficients of linear operators. 
The first linear equations are 

N 

dlS-zd[S = ^2(31^)8^, (4) 
p=i 

where d l n = -J~-, functions u^t 1 ,... ,t N ) represent first coefficients of ex- 
pansion of functions S p (z; (t 1 , . . . , t^)) at infinity; S p = YlnLi SnZ~ n , u p = 
—S p . Taking compatibility condition of a pair of linear equations @ la- 
belled by superscripts i and j, we get a closed equation for the vector field 

a£flju - did{u - [0fu, d{u] = 0, (5) 

where [u, v] is a standard commutator of vector fields. In a similar manner 
one can construct higher linear problems (j3J) and nonlinear systems (jSJ. 

Lax pairs of the type (jlj and corresponding systems © are known for 
a long time (see e.g. |23j). The method to solve Cauchy problem and the 
dressing method for Q based on the spectral theory of operators of the type 
(jlj has been recently developed by S.V. Manakov and P.M. Santini [2"1| . 

In the simplest case N = 1 our scheme provides us with the (9-dressing, 
linear problems and nonlinear equations for the universal hierarchy of hy- 
drodynamic type studied in |251 126j . (9-dressing for this hierarchy will be 
studied elsewhere. In this paper we will concentrate on the Hamiltonian 
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reductions of the equations for vector fields and associated heavenly type 
equations. 



3 Heavenly equation 

Let us consider two-component case (N = 2) and denote x = t\, y = t\, 
is easy to check that the system admits the reduction 
to Hamiltonian vector field u = Q y d x — @ x d y , for which it becomes a second 
heavenly equation 

Q t y ~ Q ix - Ql y + e xX Oyy = 0. (6) 

Hamiltonian reduction of the system © surprisingly corresponds to the 
Hamiltonian reduction of the <9-system Characterization of this re- 

duction in terms of problem Q is given by the following statements (for 
simplicity we formulate them for the two-component case, TV-component 
case is completely analogous). 

Proposition 1 Let the d-data for the problem 0) be of the form 

W 1 = W :2 , W 2 = -W t i, (7) 

where W(z, z; S , S 2 ) is some function (potential for the d-data). Then the 
two-forms 

w = 5S 1 A5S 2 , (8) 

or, equivalently, 

w{5S, SS) := SS^S 2 - SS^S 2 , (9) 

where 5 and 8 denote arbitrary variations, are analytic inside the unit disc 
(in general case in C \ G, which is the support of d-data). 

Proof. The problem (Q) in this case looks like Hamilton equations with 
complex time z and Hamiltonian W, 

- i dW - o dW 

ds = ds^ ds = ~W' (10) 

and in complete analogy with standard Hamiltonian mechanics, we get (ana- 
logue of Liouville theorem, see e.g. [2*7] ^ 

d(5S 1 A 5S 2 ) = 0. (11) 
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It is also a simple check that 

BiSS^S 2 - SS l SS 2 ) = 0. (12) 

□ 

Analyticity of the two-form w (jUJ) inside the unit disc readily implies that 
vector field S\d x + S 2 d y is Hamiltonian. 

Proposition 2 Let the d -problem 0) be of the Hamiltonian form il(J\) . 
Then the vector field S\d x + S\d y is Hamiltonian. 

Proof Identity means that for any pair of solutions f , f of linear equa- 
tions © inside the unit disc one has 

<9(/7 2 - /V 2 ) = 0. (13) 

Taking a pair of solutions f = (S x , S 2 ), f = (Sy , S 2 ) and considering analytic 
properties of the form Q, we come to the conclusion that 

{SijS^} := S\S 2 — SyS 2 = 1. (14) 

The fist term of expansion of this relation at z — > oo gives the identity 
d x S{ + d y S 2 = 0. Thus 5| = — 6^, Sf = Q x , meaning that the vector field 
S%d x + Sfdy is indeed Hamiltonian. □ 

Identity <|13|) can be used in many different ways. We will show now 
that it directly leads to the linear problems of the type PJ. Indeed, let us 
chose f = S x and f = CS, where £ is a first order differential operator in 
independent variables depending on z such that (CS)(z) is bounded outside 
the unit disc and (CS)(z) — > as z — > 0. Then due to ((T3|) w(S x ,£S) 
is analytic in the whole complex plane and vanishes at z — > oo. Hence 
w(S x , CS) = and consequently CS — aS x = 0, where a is some function 
of independent variables and z. On the other hand, choosing f = S y and 
f = £S, we come to the conclusion that CS — (3S X = 0. Then it is easy to 
see that a = j3 = and CS = 0. So we have obtained linear problem for S. 

For the heavenly equation © one gets the known linear problems (see 

e.g. m) 

- zS x - {S,@ x } = 0, . , 

S t --zS,-{S,G J/ } = 0. 1 0) 

Let us emphasize that our derivation of linear problems in the Hamil- 
tonian case (jlOjl was done for arbitrary 9-data (the only thing we need is 
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existence of solution), and we didn't suggest validity of an analogue of Liou- 
ville theorem for the system (j3J) corresponding to Hamiltonian case (in fact 
in the process of derivation we have proved it). 

4 r-function and addition formulae for the heav- 
enly equation hierarchy 

Analyticity of the two- form w (JSJ) inside the unit disc or equation (|13|) play a 
role of Hirota identity for heavenly equation hierarchy. This property can be 
formulated in a standard way for the boundary value of w on the unit circle, 
using a projection operator. It implies the existence of the r-function for the 
heavenly equation hierarchy (which coincides with introduced above for 
the heavenly equation) and gives 'addition formulae' (generating equations) 
of the hierarchy and heavenly equation itself. 

Proposition 3 Identity ill)) implies that one-form 

° = 1 hi j>(& SS l ~ ^SSfodz (16) 

is closed. 

The proof is by simple direct calculation. 

We define a r-function 0(ti,t2) for heavenly equation hierarchy through 
closed one- form ()16|) by the relation 50 = 9. Introducing vertex operators 
^(z) = J2^Li z ~ n dn, D 2 {z) = J2^=i z ~ n dni ^ is easy to demonstrate that 

S\z) = -D 2 (z)Q, S 2 {z) = D l {z)Q. (17) 
Substituting this representation into (|14[l. we get the equation 

D 2 (z)Q x - D 1 (z)Q y - {D 1 (z)0, D 2 {z)Q} = (18) 

The first nontrivial order of expansion of this equation at z — > oo gives 
exactly the heavenly equation ©. 

To derive addition formulae (generating equations in terms of vertex 
operators) for 0, we will consider the two-forms w(D 1 (z')S(z), D 1 (z")S(z)), 
w(D 2 (z')S(z), D 2 (z")S(z)), w(D 1 (z')S(z), D 2 (z")S(z)). Taking into account the 
identity ()13l) and analytic properties of these forms, one gets 

D 1 {z')S 1 {z) ■ D 1 {z")S 2 {z) - D 1 {z")S 1 {z) ■ D 1 {z')S 2 {z) 

= -J—D\z")S 2 (z>) - -^D\z')S 2 {z'% 
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D 2 (z')S l (z) ■ D 2 (z")S 2 (z) - D 2 (z")S\z) ■ D 2 (z')S 2 (z) 



^—D 2 {z')~S\z") - -J—D 2 (z")S 1 (z'), (19) 



D l (z)S l (z) ■ D 2 {z")S 2 {z) - D 2 (z")S\z) ■ D\z')S 2 {z) 

= -j 1 ^— + -^—D 2 {z")S\z') + -J—D\z')S 2 (z"). 

z' — z z" — z z' — z z" — z 

These relations directly imply the existence of the T-function (that is equiv- 
alent to the definition in terms of closed variational one- form (|16|) ) and 
provide us with addition formulae for it. Indeed, taking into account the 
asymptotic behaviour of both sides of relations (|19|) (or using the analytic 
properties of the forms w in the l.h.s. of these relations and Cauchy formula 
for the domain G), one gets 

D 1 (z")S 2 {z') - D l {z')S 2 {z") = 0, 
D 2 {z")S 1 {z') - D 2 (z')S 1 (z") = 0, 
D 2 {z")S 2 {z') + D 1 (z')S 1 (z") = 0, 

that implies ((T7|) . 

Substituting representation (|T7|) to relations ffTTTf) . we obtain a set of 
addition formulae 

—^—D 1 {z"){D 1 {z') - D 1 {z))Q - —}—D 1 {z'){D 1 {z") - D 1 {z))Q 

z — z z — z 

= D l {z")D 2 {z)Q ■ D 1 (z')D 1 (z)e - D 1 {z")D 1 {z)Q ■ D l (z')D 2 {z)Q, 



-^—D 2 {z'){D 2 {z") - D 2 (z))Q - -^—D 2 {z"){D 2 {z') - D 2 {z))Q 

= D 2 {z')D 2 {z)Q ■ D 2 (z")D 1 (z)e - D 2 (z')D l (z)Q ■ D 2 (z")D 2 (z)@, (20) 



—D 2 (z")(D 1 {z') - D 1 {z))Q 1 —D 1 (z')(D 2 (z") - D 2 {z))Q 



= D\z')D 1 (z)Q ■ D 2 (z")D 2 (z)Q - D 1 {z')D 2 {z)Q ■ D 2 {z")D l {z)Q. 

Expansion of these equations into powers of parameters z, z", z" generates 
heavenly equation hierarchy. 
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5 <9-dressing method and r-function for the heav- 
enly equation hierarchy 

Similar to the case of dispersionless integrable hierarchies |18l 119] . it is pos- 
sible to obtain explicit formula for the r-function of heavenly equation hi- 
erarchy, which is given by the action for the system (jlUj) evaluated on the 
solution of this system. 

The problem 1)10(1 can be obtained by variation of the action 

f = —ff (^BS 1 -W{z,z,S\S 2 )) dz Adz, (21) 

2vriyy C \G v 1 

where one should consider independent variations of S, possessing required 
analytic properties (analytic in G, decreasing at infinity), keeping So fixed. 

Proposition 4 The function 

0(t) = ^-JJ (s 2 (t)dS l (t) -W(z,z,S\t),S 2 (t))^j dz Adz, (22) 

where t = (ti,t2), D is a unit disc, i.e., the action \21\) evaluated on the 
solution of the problem U(J\) . is a T-function of the heavenly equation hier- 
archy. 

Proof. Considering @(t) as a functional of So and calculating its variation, 
we get 

50 = — (f{S 2 5Sl - S^S^dz, 
2tti J 

that coincides exactly with one-form (|16|) used to define the r-function. □ 
Let us consider also the action 

f = —JJ (S 2 dS l -W(z,z,S 1 ,S 2 ))dz Adz, (23) 

which is an exact analogue of the classical action associated with the Hamil- 
tonian equations (fTU|) . Then 

f = ±.Jsfs^ + f. 

So on the solutions of the <9-problem (pT)|) one has (0 = /) 

e = e-VV— 

n=l " 
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Similar to dispersionless case [181 I19j . the formula (122() provides us with 
infinitesimal symmetries for the addition formulae (|2Uj) and in particular for 
the heavenly equation (JfjJ). They are given by the expression 

5Q = — — [f SW(z,z,S)dz Adz, (24) 
2vri JJ D 

where 5W(z, z, S) is an arbitrary function in the unit disc. 

In the simplest case 8W = 5(z — zq)F(S(zq, t)), where F is an arbitrary 
function, one has 

<ye = F(s(*d,t)). (25) 

It is a simple direct check that an arbitrary function F(S 1 ,S 2 ) of the solu- 
tions S 1 and S 2 of linear problems (|15|) solves the linear equation 

(SQ) ty - (S@) ix + {9,, (50) y } - {Q y , (5Q)x} = 0, 

which defines infinitesimal symmetries of the heavenly equation © . 

Finally we note that, similar to the classical mechanics, canonical trans- 
formations S l — ► S n = p(S 1 ,S 2 ) preserve the two- form w (jSJ) and leave 
9-problem ()1U|) invariant. In addition, associated linear problems CS = 
are invariant under these transformations. 

Under canonical transformation with the generating function F one has 
W — > W' = W + Fg. So for canonical transformation independent of z the 
G-function corresponding to (|23|) and r-function ()22() remain invariant while 
in general case G' = Q + ^ § Fdz. 



6 Conclusion 

The results presented above are generalizable naturally to multicomponent 
case. Corresponding systems admit various reductions. In particular, for 
even number of components (2iV), the Hamiltonian reduction, 

N 

8=1 

where Xi = tj, % = 1, . . . , 2N, gives rise to the equations 

®Xit k ~ ®x k ti ~ {Qxi,&x k } = 0, 
where X{ = t 2 , i = 1, . . . , 2N, and 

ir-ir-)- 

^ \OX2p OX2p-l OX 2p -l OX 2 pJ 
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These equations are associated with linear problems 



Stj - zS Xi = {S, @ Xi }. 
The reduction u = ^« + v ) where 

N 

p=M+l 

leads to equations 

e^-e^- {6^,9^1 = 0, (26) 

where the Poisson bracket is of the form 

p Jm+i V^sp Cx 2p -i ax 2 p~i ox 2p 

In the case iV = 2, M = 1, taking i = 1, k = 2 one obtains six-dimensional 
generalization of the second heavenly equation [5]. The (^-dressing method 
for these multidimensional equations will be studied in elsewhere. 
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